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Abstract 

We prove new identities for the constant in the Mertens product over primes in the 
arithmetic progressions a mod q. AMS Classification: 11N13 

1 Introduction 

Let a, q be integers with {q,a) = 1 and denote by p a prime number. In 1974 Williams |[8l 
proved that 

-p-f / 1 \ C{q,a) f ^ \ 

P{x;q,a)= n ll-^j = (iog.)i/cp(.)+^^(log.)iM.mJ 

p=a mod q 

as — s> +00, where C{q, a) is real and positive and satisfies 



cp(^)/^i„VL(l,5c)^ 

where y is the Euler constant, cp is the Euler totient function, L{s, %) is the Dirichlet L- function 
associated to the Dirichlet character % mod q and %o is the principal character to the modulus 
q. The function K is defined by means of 

-|-oo 

n=l 

where k-^{n) is the completely multiplicative function whose value at primes is given by 

In our recent paper Pl we obtained a version of Williams's result stated in ^ which is uniform 
in the q aspect. In the same paper, as a by-product, we also obtained the following elementary 
expression for C{q, a) : 

c(,,.)*<")=.-vn(i-i)*"°' (2) 
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where a{p;q,a) =(p{q) — I if p = a mod q and a{p; q,a) = —1 otherwise. The infinite product 
is convergent, though not absolutely, by the Prime Number Theorem for Arithmetic Progres- 
sions. Actually, a slightly simpler proof of identity ^ than the one we gave in Q, sect. 6, can 
be obtained as follows, once the relevant limit is known to exist: taking logarithms in ([T]) and 
using the classical Mertens Theorem, we find 



logC(^,a)=^lim {log n 



p=a mod q 



p=a mod q 

— ^- + — -log<^ lim n 1-- 



The product in (|2l) is very slowly convergent and it is difficult to compute an accurate numerical 
approximation to C(^, a) from it. Our aim here to give a different form for the constant defined 
in ©I unfortunately, this form is not suitable for numerical computations, a problem we tackle 
in part II [5]. 

As a corollary of the identities proved in the first part of the paper, we derive the formulae 
that Uchiyama [|7l gave in the case q = A and a E {1, 3}, though Uchiyama's direct proof is 
obviously much simpler. We also derive the explicit expressions that Williams gave in Theo- 
rem 2 of [jSJ for C(24, a) for every integer a such that (24, a) = 1 and the ones that Grosswald 
[[3| obtained for C{q,a) for q E {4, 6, 8} and every integer a coprime to q. We recall that, in 
Proposition 1 of flU, Moree gives formula ^ below when ^ is a prime number and a= I. 

The statement of these new formulae themselves is not simple, and reflects both the struc- 
ture of the group Z* and the properties of the residue class a. For this reason, we will not state 
a formal theorem here, but rather point to the various results as formulae ([8]) or its alternative 
version Q when Z* is cyclic and a = I, then (fT3l) for general q and a = I, and (fT5l) in the 
most general case. We think our results will be clearer if we go through stages of increasing 
generality. 

We may summarize our main result saying that, in the case a= 1, for each reduced residue 
class b we have to determine a positive integer tt (actually, the order of b in the multiplicative 
group Z*) and we will then express C{q, 1) as a sort of Euler product where a prime p has the 
exponent —tp. Collecting all residue classes of maximal order, we may reduce the number of 
factors needed, at the price of the computation of the power of a suitable value of the Riemann 
zeta function at an even integer. The important feature of (fT3l) is that t^, >2 for all b 1, 
whereas the exponent of all primes in the Euler factors in (|2]) are — 1 . 

The case I is genuinely more complicated: in fact, in general it will not be possible to 
give a simple closed form for the Euler factors, though they can be expressed by means of a 
rapidly convergent power series. On the other hand, the constant defined in (fT6l) will arise: it 
is related to the so-called Meissel-Mertens constant; for its computation in the case q = I see 
§2.2 of Finch [2J. 

We would like to thank Pieter Moree for providing us some references and Giuseppe 
Molteni for suggesting a simplification in the proof of Lemma [H 
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2 Reduction to character sums 



It turns out to be better to get rid of the prime factors of q at the outset: therefore, we let c{q, a) 
be defined by means of 

C(^,a)W) = g-Y_^c(^,a) so that c[ci,d) =\\\\ - -\ (3) 



1 \ 



p\<l 

and let 



c{x;q,a) = Y[{^ 



p<x ^ P 

denote its partial product. Our strategy is to express the product of c{x; q, a) and partial products 
of powers of L(l,x) as a complicated but quickly convergent product, where % ranges over all 
non-principal Dirichlet characters modulo q. When necessary, we use the abbreviation 

U{q,a)= n (4) 

X mod q 

Xt^Xo 

By orthogonality, we have 

X mod q p<x ^ p<x x mod q ^ ^ 

Xt^Xo p\q xt^Xo 



Using the Taylor series expansion of log(l — t) we see that 



Xmodg ^ P ^ ^ m>2'^P 



(6) 



Xt^Xo 

say, where Smip'^q^o) is the character sum defined by 

S,n{p;q.a)= m{x"'{p)-l{p)). (7) 

X mod q 

It is this character sum that reflects the structure of Z* and the properties of the element a. 
We will prove below in (flOl) and (fTTj) that either p = a mod q, or Smip^q^ci) vanishes unless m 
belongs to a suitable arithmetic progression modulo a divisor of cp(^). The simplest case, not 
surprisingly, is when Z* is cyclic and a = I. 



3 The character sum S,n in the simplest case 

We notice that, obviously, S,n{l;q,a) = 0, and we may assume that p^ \ mod q. For the time 
being, we also assume that a = I. Let tp denote the order of p in the multiplicative group 
Z*, that is, the smallest positive integer k such that p^ = I mod q, and notice that tp>2 since 
p ^ 1 mod q. It is then quite easy to see that 



I %'"{p)= I xip'n 

X mod q X mod 1 




if tp\ m 

Otherwise. 
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Hence, using orthogonality and the Taylor series for log(l —t) again, we have 



m>2 



-Sm{p;q, 1 = L , nt„ = 1 ~ ^ 

mp™ ntpp p \ p i' 



We classify primes according to their residue class b mod q, and notice that tp depends only on 
b,if p = b mod q. Substituting into ^ and letting x +o°, we see that 

c(.,i)^n(„,,) n n (8, 

p=b mod q ^ 

We notice that the quantity Il{q,l) is connected to the Dedekind zeta function of the q-th 
cyclotomic field ^ = Q(^^) by means of the relation 

Xf{p)\ 



n(^,iri=Resc^(.) n n(i-^) 



X mod q p\q 

where Xf denotes the primitive character that induces % and / is its conductor. 

Assume that Z* is cyclic. Relation ([8]) is our first formula, and it is worth noticing that a 
slightly better form, from the point of view of the explicit computation of C{q, 1), can be given 
grouping the contribution of primes of maximal order. More specifically, we may rewrite ([8]) as 



c(,,i)=n(,,i) n n(i-^ 



nW(q) p 

n>\ tp=n 



-V?{q)/n 



p ^ p ^ 

tp=(?{q) l<n«?{q) tp-" 

X n{(.-^)-"'""(.-^)} 

nmq) p 
l<n<(p(g) tp~" 

c(cp(^))n(,,i)n(i--^) _n (i 



p'^ 



X 



p\q p=l mod q r 

n, n{(.-^)-^"""('-^4T)}- 

"l<p(?) p ^ ^ 

\<n<<S){q) tp-" 

The value of C,{(p{q)) is easily computed, at least when q is comparatively small, by means of 
the Bernoulli numbers, since (p{q) is even for q > 3. It is also worth noticing that the expo- 
nents of the prime numbers in the last product are all at least 2, though they are usually much 
larger. Uchiyama's formula (fTTI) for C(4, 1) in |I71 is the case ^ = 4 of the above expression: 
there is only one non-principal character % modulo 4, and L(l,%) = 7l/4. The last product is 
empty, and ^(2) = 71^/6. The formula for C(4,3) is easily deduced from this and the classical 
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Mertens Theorem since C(4, 3)C(4, 1) = 2e^^. Grosswald's fomiula (fTSi) for C(6, 1) in [3] is 
another special case of © and Q since (p(6) = 2, there is only one non-principal character % 
modulo 6, and = 7r/(2v^). By the Mertens Theorem we have C(6,5)C(6, 1) = 3e^^ 

and the formula for C(6,5) can be easily deduced. Our formula dH) also contains Moree's |l6l 
which is the special case where q is prime. 



4 The sum Sm in the general case 

Lemma 1 Let Sm be the character sum defined in dV]). If a = b mod q then there exists a positive 
integer ta dividing cp(^) such that 



S,n[a;q,a) = { . (10) 

U otherwise. 



In this case, ta is precisely the order of a in the group Z*, so that ta >2 unless a = b = 1. 
If a ^ b mod q then either the equation b^ = a mod q has no solution y ^N, or there exist a 
positive integer t}, dividing (p{q) and an integer s^ such that i < st, <th and 

I (p(q) ifb^' = a mod q has a solution, and m = Sb mod tt,, 
Sm(b;q,a) = < . (11) 

I U otherwise. 

In particular, if a = 1 then 5m(l;<?, 1) = 0, while, ifb ^ 1 mod q, then there is an integer th>2 
such that Sm{b;q, 1) = (p{q) ifm = mod t}, and is otherwise. 

Proof. The complete multiplicativity of the Dirichlet characters implies that 

Sm{b;q,a)= £ {x{b'"a-')-x{ba-')) 

X mod q 

If a = b = I mod q then S,„{ci',q, a) = 0. If a = b ^ I mod q then the equation a^ = a mod q 
has the solution y = I mod ta, where ta is the order of a in Z*. Hence, S,n{a;q,a) = if 
m = 1 mod ta and S,n{a',q-,ci) = — cp(^) otherwise, by orthogonality. 

If a ^ b mod q and the equation b^' = a mod q has no solution, then Sm{b,q,a) = by 
orthogonality. If the equation above has the solution y = sy mod ty (where > 2 is a suitable 
divisor of (p((?) and Sh is an integer with 1 < < tb) then S,n{b;q,a) = (p{q) if m = Sb mod tb 
and is otherwise. □ 

It is quite clear from the Lemma above that the case a = 1 is indeed much simpler than the 
general one. We see that we have proved that ® holds also in the general case, where tb is the 
divisor of ^){q) such that the solution of the equation M = 1 mod q is the class mod tb, that 
is, the order of b inZ*. The computations that lead to ^ are still essentially valid, with one 
modification. We recall the definition of the Carmichael X function: if ^ = 2"p"' ■ ■ • p"*, then 

ricm{(p(2«),(p(;.«>),...,(p(;.«^)} ifa<2; 
^^^^ \lcm{2«-2,cp(;.?"),...,cp(;.«^)} ifa>3. ^ ^ 

In other words, ^{q) is the highest order of the elements of the group Z*. Let .^(^) = {b E 
Z* \ {1} : b^ = 1 mod q for some positive k < ^{q)} denote the set of elements of Z* \ {1} 
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whose order is not maximal. Arguing as in the proof of dH), that is, grouping the contribution of 
the primes of maximal order, we obtain the following Theorem, which generalizes Williams's 
ifSl and Grosswald's jSl formulae for C(24, 1) and C(8, 1) respectively and contains Q as a 
special case. 

Theorem 1 For all integers q>3 the value of the constant c{q, 1) is given by 

c(„ .) ^ n(„ 1, n(. - ^r'^'" n (> - i)"""" 

p\q " p=\ mod q " 

n,__n,{(i-^) (i-^) } (13) 



X 

b<EJ^ {q) p=b mod q 



where n(^, a) is defined in dU, X is the Carmichael lambda function defined in (fT2)) . t}, denotes 
the order ofb in the group Z* and ^ (^) = {Z? G Z* \ { 1 } : t^ < X^q) }. 

5 The general formula 

In this section we assume that a^ \ mod q, and we let ® {q, a) denote the set {Z? G Z* \ { 1 , a} : 
the equation M = a mod q has a solution}. For the elements of this set, we implicitly define 
the integers t^ and sy as in the proof of Lemma [IJ In order to state the formula corresponding 
to (fT3l) in general, we need to introduce the function ft,s which is defined, for positive integers 
t and s with \ <s <t and real x with |x| < 1, by means of the relation 

n=s mod / 

The rightmost equality is proved computing the derivative of the function f,s and then summing 
the ensuing geometric progression. Notice that, when s = t, di closed form for the integral 
can be easily given in terms of the logarithmic function, as we did above: indeed, ft,t{x) = 
—t^^ log(l —x^). In general, the Taylor series for f^s is fairly quickly convergent since we will 
compute it at X = . 

Lemma [T] above amounts to saying that, given integers q and a and a reduced residue class 
b mod cp(^), either the equation b^ = a mod q does not have a solution, or its solution is a 
congruence class Sb mod tt,, where tf, \ (p(^) and we may assume that 1 < ^z, < and that tb > 2, 
since th = I if and only if a = b = 1. Moreover, unless a = b, the congruence class Sb will not 
be 1 mod q. Therefore, classifying primes according to their residue class b mod q again, and 
substituting either (fTOl) or (fTTI) into we see that the corresponding factor in the product © 
is I if b = I or b ^ {q^a) U {a}, and is 

-1- 



n exp(^cp(^)4,,,(;? 

p=b mod q 

ifbE'B {q, a) . Recall that for b = a mod q we have Sa = i: hence, for primes p = a mod q we 
have a factor 

n exp(-(p(^) I, U exp((p(^) £ 

p=a mod q ^ m>2 ^ p=a mod ^ ^ m>2 '"^^ ra>2 '^"^^ ^ 

1 mod ta m= 1 mod ta 
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Finally, collecting all identities, we see that we have proved the following result. 

Theorem 2 For all integers q>3 and all integers a such that {q, a) = \ and a ^ 1 mod q, the 
value of the constant c{q, a) is given by 



n(^,a) is defined in dH) and ft^s is defined in (fT4l) . 

For the special cases q G {4, 6, 8, 24} and {q,a) = 1 with a ^ I mod q, equation (fTSi) 
collapses to the formulae given by Uchiyama [7J, Williams dH and Grosswald 

For some special values of t and s it is possible to compute a closed form for ,,, as in 
the previous sections, and give a more explicit result: the following section contains some 
examples. 

6 Explicit values 

Using ([3]), ^ and (fT3l) . we can compute explicitly a few values of the constant C{q, a) . For the 
evaluation of L(l,x) needed to determine n(^, 1) we refer to Corollary 10.3.2 and Proposition 
10.3.5 of Cohen The value (flTI) is due to Uchiyama while the values (fTSi) and (fT9l) are due 
to Grosswald. Notice that Jl{S) = and that n(8, 1) = 327r"2Qog(3 _^ 2V2))-\ 




p=a mod q 



(15) 



where ?B a) = {Z? G Z* \ { 1 , a} : the equation by = a mod q has a solution}, 




(16) 




p= 1 mod 4 



(17) 




(18) 



C(8,l)4 



32 



1 



71^^-^(8,1) _n (i 



/:)=1 mod 8 




log(3 + 2y2) 



p= 1 mod 8 
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For ^ = 24 we recover the value given on page 357 of Williams [[HI: using ([3]) and ([T3l) we get 

c(24i)«= n fi--'' 

91og(2 + v^)log(l + v^)log(5 + 2V6)p,iUd24^ 
since (p(24) = 8, X(24) = 2, (24) = and 

n(24,i) 



7r4log(2 + v^) log( 1 + V2) log(5 + 2^6) ' 
For the more complicated case q= 15 we get 

p=b mod 15 

where we have used the fact that J? (15) = {4, 11, 14} and the values ^(4) = 7r'^/90 andX(15) 
4. A fairly lengthy computation reveals that 

,x 1 2^71^^ / 9/ sin(7r/15) \ , 9 /sin(77r/15)\\ , /l + VS 

This can be made more explicit using suitable trigonometrical identities and the value for 
sin(7r/15), 

Using ^ we can compute C(5, 1): 

c(v)^^i^.-n(5,,)|i n (.-^) n ('-^)"('-^^ 

^ "^-^p^lmodS " ^ p=4mod5 " 

150 log((l + v^)/2) 

mod 5 P p=4 mod 5 ^ ^ 

where we have used the value n(5, 1) = 25-\/57r^2(^4iog((l + -\/5)/2))^^ and the fact that 
^(5) = {4}. 

The last examples are for ^ = 5 and a G {2, 3, 4}. A short computation shows that 



/4,i(x) = -log(^Y3^j + 2^^^^"W 

/4,2(x) = ^l0g(^- 



-x2 



^ , , 1 /I +X\ 1 , , 

Furthermore, using the fact that is generated by 2, we see that (B (5, 2) = {3} (with ^3 = 3 
and ?3 = 4), S (5, 3) = {2} (with ^2 = 3 and t2 = 4) and -B (5, 4) = {2, 3} (with ^2 = ^3 = 2 and 
f2 = ?3 = 4). These results show that 

c(5,2)=n(5,2)exp(45(5,2)) H exp(4/4,3 (-)) exp(4 £ UJ^-l. 

p=3mod5 \p// \ p=2mod5 P P 
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p=b mod 5 



and the corresponding values for the constants C(5,a) can be found using ([3]). In the case of 
C(5,4) we can be slightly more explicit since/2, i(x) = l/21og((H-x)/(l — x)), so that 

c(5,4)=n(5,4)exp(45(5,4)) Y[ {\^) Y[ (|^)V^/^. 

p=b mod 5 
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